function grad_b_value = compute_grad_b()
% b1 = @(x,y) 1 - 2.*x;
% b2 = @(x,y) 2.*x + 2.*y -1;
% b3 = @(x,y) 1 - 2.*y;
b1_x = @(x,y) -2 + 0.*x;
b1_y = @(x,y) 0.*x + 0.*y;
b2_x = @(x,y) 2 + 0.*y;
b2_y = @(x,y) 2 + 0.*x;
b3_x = @(x,y) 0.*x;
b3_y = @(x,y) -2 + 0.*y;
grad_b_value = zeros(3,3,2);
gauss_points = [0.1666666666667, 0.1666666666667;
  0.6666666666667,0.1666666666667;
   0.1666666666667,0.6666666666667];
for i = 1 : 3
    grad_b_value(1,i,1) = -2;
    grad_b_value(1,i,2) = 0;
    grad_b_value(2,i,1) = 2;
    grad_b_value(2,i,2) = 2;
    grad_b_value(3,i,1) = 0;
    grad_b_value(3,i,2) = -2;
end
end